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线性 乘 性 规划 的 全 局 优化 算法 
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(1- 中 南大 学 数学 科学 与 计算 技术 学 院 ， 长 沙 410075; 2- 广东 金融 学 院 应 用 数学 系 ， 广 州 510521) 

摘 要: 本 文 研究 线性 乘 性 规划 问题 (LMP) 的 全 局 最 优化 算法 ， 线 性 乘 性 规划 问题 在 生产 运输 、 工 厂 布 
局 设计 、 超 大 规模 集成 电路 芯片 设计 等 方面 有 重要 的 应 用 。 首 先 将 LMP 问题 转化 为 等 价 规划 问 
题 (P1)， 然 后 利用 参数 线性 化 方法 在 相应 的 超 和 矩形 上 求 得 问题 (P1) 的 目标 函数 和 约束 函数 线性 
下 界 估计 ， 并 提出 了 一 个 求 线性 乘 性 规划 全 局 解 的 确定 性 全 局 优化 算法 ， 并 证 明了 算法 的 收敛 
性 。 数 值 实验 表明 提出 的 方法 是 可 行 和 有 效 的 。 
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1 引言 


考虑 下 面 线性 乘 性 规划 的 全 局 优化 问题 
min 2Z=Go(z) 
LMP(X?) st. G;(£) <ð; j-.12.-,m, 


z€X?-[r,z|c RN, 
其 中 
Tj N 
Gj) = [| ( 32s + di), j=0,1, ,m, 
i-1 n=l 
Yitn, djt 都 是 任意 的 实数 
N 
oa mscda20, $20, j201--,m t-12.-,T, z»0, g«z 
Hel 
都 是 有 限 实 向 量 。 

上 述 线性 乘 性 规划 广泛 应 用 于 生产 运输 踢 、 工 厂 布局 设计 四、 超大 规模 集成 电路 芯片 设 
计 回 、 数 据 挖掘 与 模式 识别 向 等 实际 问题 。 尽 管 求解 LMP 的 优化 方法 很 多 ， 而 对 LMP 的 全 局 
优化 方法 却 很 少 ， 已 有 的 全 局 优化 方法 都 是 针对 LMP 的 特殊 形式 进行 讨论 的 上 F999。 本 文 给 出 
一 种 参数 线性 化 枝 术 ， 对 LMP 问题 提出 了 一 个 分 支 定 界 全 局 优化 算法 ， 分 支 定 界 算法 是 一 种 
比较 有 效 的 全 局 优化 算法 ， 在 很 多 全 局 优化 问题 都 有 重要 的 应 用 [i913 。 利 用 对 数 变 换 及 目标 
函数 和 约束 函数 的 线性 下 界 估 计 ， 建 立 LMP 的 线性 松 驰 规划 问题 LRP， 通 过 求解 一 系列 线性 
松 驰 规划 问题 LRP， 逐 步 改进 原 LMP 最 优 目 标 值 的 上 界 和 下 界 ， 最 终 确 定 原 问 题 的 全 局 最 小 
解 。 

收 稿 日 期 : 2008-07-18. 作者 简介 : 周 雪 刚 (1974 年 1 月 生 )， 男 ， 博 士 ， 讲 师 . 研究 方向 : 全 局 优化 与 模糊 优化 . 
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2 ”线性 化 技术 


为 求解 LMP， 对 目标 函数 与 乘 性 约束 不 等 式 取 对 数 ， 则 LMP 可 转化 为 如 下 等 价 的 规划 问 
题 








min WVo(z) 
P(X?) st. V;(z)&lnó; j-12.-,m, (1) 


x € X? = |z?, Z°] c RN, 
其 中 对 所 有 的 j= 0,1, ,m， 有 
Tj N 
Vi(z) = > in( buon Æ dje). 
t-1 n-l 
设 和 =[fz, 姑 ={zeRNIz 和 zz 入 下， 对 任意 的 ze 和， 定义 
z(c)-z-o(s- mz), (2) 
其 中 oc € {0,1}N 是 NN 维 向 量 ， 其 分 量 o; = 1 或 0。 显 然 ，z(0) = z, z(1) = z, 对 任意 
Ho € {0,1}v，z(c) 表 示 区 间 向 量 瑟 的 一 个 顶点 。 下 面 定 理 说 明 如 何 构造 函数 亚 j(z) (7 = 


0,1,.… m 的 上 、 下 界线 性 函数 。 
定理 2.1 对 于 任意 的 区 闻 向 量 XC X? c RN HEX EKAS 


Yi(z) = ln ( S T de) 
i-l 


及 其 梯度 





voo = (8:82... 4497 


则 存在 z, Z € RN 及 线性 函数 


VWi(z; X,o) := z(o)” -x + Vi(a(o)) - z(o)” - a(o), (3) 
V" (z; X,0) := zl — 0o)" -x + Wi(z(o)) -all -—0o) : a(c), (4) 
对 所 有 的 ze 和 满足 
Wi(z; X,o) < W(z) < Vt; X, o), (5) 
V (z(o); X, v) = V? (z(o); X, c) = Y; (z(0)), (6) 


其 中 z(o)，z(c) 分 别 表示 区 间 向 量 X 和 及 := [ze Z 02N AUS al oa) PHE- o= 
(1 一 01 , 工 一 IN)7。 
证 明 X DECIR] Fd EE X HW (a), FENE z = (za szin)! M Z = (u.c Zen)? 满 
A 
z S V(r) SU, VIEX, (7) 
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其 中 
N ed N =i 
Zt = min fa 人 5 min (37tnZn; imn) + 4) , | 5y max (^m, Ynzn) + de) }, 
n-l n=1 
N _1 N ze 
Zi = max EA 》 min (Ynzn, YnEn) + de) i 人 ( 》 max (Ynzn, inn) 十 d) y 
n=l n=1 


i 二 1,2,.… N, m 表示 zz 的 第 i 个 分 量 。 
由 微分 中 值 定理 ， 对 任意 zeE 和 有 


V,(z) = V.(z(o)) + YE)T - (z — z(0)), 
其 中 = jz 十 (1 一 pj)z(o), ef0,1。 根 据 式 (7)， 对 任意 的 ze 和 ， 若 ui = 0， 则 有 不 等 式 





cm > zulo), zi— zi(0) > 0 
PA < zulo), m-m()«0 
成 立 ， 所 以 不 论 是 mi = 0 还 是 0 = 1 都 有 
PUO (a — 19) > zulo) (es = nin) 


成 立 。 因 而 有 


Vi(z) = )* > Mm - zi(0)) 


WV 


N 
Vr(z(o)) + Y zulo) (zi — zi(o)) 
i—1l 


= alo)" -x+ Vi(z(o)) — zlo)" - z(o). 


所 以 对 任意 的 ze X, A Viz X, o) < V2) M Vi(z(o); X, o) = Wi(z(o)) RE. 
同 理 ， 对 o; = 0， 则 有 不 等 式 
89, 9v. (£) €zu(l—o0), zi—z(co)20, vrexX 
(O6 
RL, HHP zll- a) 表示 zz(1 一 o) 的 第 i 个 分 量 。 若 o; — 1， 则 有 不 等 式 


2z(1-o) Zi~zilo) <0, vreX, 





所 以 不 论 是 o; = 0 还 是 o; = 1 都 有 


OW:(é) 
a&i 





(zi — zi(o)) < zu(1 — o) (zi — zi(o)) 
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成 立 。 所 以 由 中 值 定理 可 得 


V,(z) = ) 十 Se e UMS. U zi(a)) 


N 
V,(x(o)) + 3 ua — e) (zi — zi(o)) 


i=1 


= zu(1— o)" -x+ Vi(z(o)) 一 2(1 — 0o)" - x(e). 


IN 


因此 ， 对 任意 的 ze X, d Wr(z; X,o) 2 V,(x) 和 更 8(z(c);Xa) = Vi(z(o)) RE. 


根据 定理 2.1， 对 任意 多 = [nz CX c RY, j=0,1, ,mm, t= 1,… Tp SAEZ 
数 o € {0,1jw， 能 按照 如 下 方式 构造 函数 更 ; : X? 一 RRETAN VF (z; X, 0) 和 线性 上 


RAA VY (z; X,c)。 对 函数 
Yj(7) = m( Dam +d) j=0,1, ,m, 
及 任 一 区 间 向 量 基 ， 计 算 Zj := [zje Zi] 使 其 满足 不 等 式 (7)， 其 中 


Za—(nacc dux Zt = (Zinc Zan), 


对 i = 1,2, N, j—0,1,-,m, t=1,.… T}, 设 


N 
Xjt = DY Yitnzn + dj, 


n-l 
并 计算 
N N 
X= » min(4;nZ,,Yinfa)tdj;, Xi = Xo Imax(TjtnZny YitnTn) + djt; 
n=} n=l 
Zjti = min nis ( XL) Xh Zji = max {a X), matX) 1) 
还 有 


N N 
z(c) - a Yon(Bn — En)en, zilo) = zji + Y. n tn — Ej )en 


n=1 n=1 


分 别 表示 区 间 向 量 X 和 221 的 顶点 。 再 根据 定理 2.1， 函 数 
N 
V jx) = In ( 20m 十 dje) 
n-—i 


的 线性 上 、 下 界 函 数 分 别 是 
VC (r; X,o) := za(1—0)7 -x + Vj(z(o)) — z«(1— c)? - z(e), 


VL (z; X,o) := zu(o)) -x£ + Vi(z(o)) - zilo)” - x(o). 


第 4 期 ASNS: 线性 乘 性 规划 的 全 局 优化 算法 647 


定义 
Tj Tj 
VL(z, X0) = 9 Yil: X 0) W (z; X,0) = ) , Vir X0). 
t=1 t=1 
显然 有 


VL (z; X, o) < Vj(z) < VU (z; X, c). (8) 

定理 2.2 B VL (n; X, o), Vc) RII Gs X, o) WERE Xegoc(0,1)*. WER X = 

[z à] c x°, rE X 和 5 = 0,1 m» ju AJ 二 VU (z; X,0) — yj) MA; 二 V;(z)— 
VL (az; X, o) 满足 


lin max A = lim max A — 0. (9) 
llz-z|—0 zEX ls—ZI—0 zEX 


证 明 对 任意 XC Xo 和 ze 大， 根据 定理 21 和 式 (1),(3), (5 与 (8) 有 


Tj N Tj 
V;(z)- n ( Y Titmzn 十 ds) = »»MPO (10) 
t-1 n=l t=1 
Tj 
VL (z; X, 0) = > 9 X0). (11) 
t=1 


根据 定理 2.1， 对 任意 的 2E XU 
0< Vj) di 9L (a; X, a) € (v5,(£) = zit(0)) (z p zf(c))， V jt (12) 


其 中 上 = jz 二 (1 pele), pE [9,1], z(a), z(o) 4 8E X M Zj, :一 [zj £j] 的 项 点 。 根据 
式 (2) 和 定理 2.1， 对 任意 固定 的 c， 不 等 式 (12) 满足 


(V^) -zxz(o)) (一 z(o)) 


N 5 
< > pm z zjmn(o)| |zn 一 zn(o)| € 3 |zji« (0) 一 zin (1|: [zn 一 rn(o)| 


n=} 


N N 
> > |Zjtn = Zitn| [zn i zn(o)| 去 5 Ital: i7 一 (x4 i |En 一 zn(a]| 
一 n=1 


n 


mm 


z 


N 

< Y henl (En — 2) (0| NX = l < (9I ( 3c had: En- zj)). 
n=l n=l 

Jtt np = pXI (0 — XAR p € [01]. 8 SER 第 二 个 不 等 式 与 第 一 个 等 式 成 


立 是 根据 式 (2) 和 (7)， 第 二 个 等 式 成 立 是 因为 对 任意 的 n= 1,2, … ,N 
lj — Zitnl = [min (ris (XI), yin COR 1 — max Cris XE) G0 H 
= meX27 eO) |= hell X307 E 
ee (2) 以 及 对 函数 z-! 关于 区 间 [X1,, X5,] 应 用 中 值 定理 ， 第 四 个 不 等 式 是 


N N 
X5, — XE = Yo [minima tin) 一 max(Yytnzn， *jtnEa)| = 》 linll lEn ~ Zal- 


n=1 n=l 
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所 以 有 
lim ax A? 一 0. 
lz- Dh,0 eX 
同 理 ， 可 以 证 明 


lim max A1 一 0. 
llz-z| «eX 


根据 定理 2.2， 当 ||F 一 zj ^ ORT, RA Vis; X, 0) 8I VU (z; X, 0) 能 充分 逼近 函数 更 j(z)。 
这 样 我 们 可 得 LMP 的 线性 松 驰 规划 问题 (LRP) 为 


min Wi(z;X,c) 
LRP(X) st. WL(zX,)«lnó;, j=1,2,...,m, (13) 
reXcx? 


设 问题 (P) 的 最 优 目 标 值 用 v(P) 表 示 ， 则 根据 上 面 的 讨论 ， 对 任意 的 和 C XO, H 
题 LMP (X) M LRP (X) 的 最 优 目标 值 之 间 有 不 等 式 v(LRP(X)) < v(LMP(X)) 成 立 。 


3 ”全 局 优化 算法 


下 面 给 出 求解 LMP 的 分 支 定 界 全 局 优化 算法 ， 算 法 的 核心 是 分 支 与 定 上 、 下 界 。 

分 支 定 界 方法 是 根据 分 支 规 则 细 分 超 和 矩形 XX0 为 有 限 个 子 超 矩 形 ， 每 个 子 超 矩 形 对 应 分 支 定 
界 树 的 一 个 节点 及 其 上 的 线性 松 驰 子 问 题 。 在 算法 的 第 k(k > 1) 次 迭代 中 ， 设 Bx 表示 由 活动 
节点 ( 即 可 能 存在 全 局 最 优 解 的 子 超 矩形 X) 构成 的 集合 。 首 先 利用 一 个 扫描 程序 将 不 可 能 包 
含 全 局 最 优 解 的 子 超 矩形 从 Qu 中 删除 ， 设 Quia 是 由 剩 下 的 子 超 矩 形 组 成 。 若 Qi+i z 0, XX 
取 Xr e Qui 进一步 考虑 。 然 后 根据 分 支 规则 ， 将 Xx 细 分 成 两 个 子 超 矩 形 Xe+1，XA+1， 完 
成 此 细 分 的 分 支 规则 称 之 为 对 分 法 。 设 X* = [rz] C X0 C RN， 对 分 法 的 步骤 如 下 : 

(i) BT -r= max(z; — Zi}, T= {1,2,.…,N)}; 


^ M . —-— * Lz 
(i) Bay WE min{z} — zj, zj- s3} = 2(2; — 2j); 
(ii) Xx 被 细 分 为 两 个 N 维 子 超 矩 形 XH, xb, 
AP = {r e R7 |z; < Ti < Ti, ij X Tj <S rj, i—-12.- N, izj) 
X$+ = {x € R7 |z; < ti < Ti, z} S Tj S Tj, i= 1,2,- ,N, ix*j). 


IERT, Qiri = Qa MXT)UUGH, X2). 

算法 的 第 二 个 基本 程序 是 定 下 界 程序 。 首 先 ，k(k > 0) 次 迭代 中 ， 对 由 分 支 程序 产生 
的 每 个 子 超 矩形 X C X0， 通 过 计算 瑟 上 的 线性 松 驰 规 划 LRP (X) 求 得 其 最 优 值 LB (X). 
iti LB (X) 即 为 问题 P1(X) 的 最 优 目 标 值 v(X) 的 一 个 下 界 。 然 后 ， 在 第 k(k > RARP, 
原 问题 的 下 界 LBx X LB = min(LB(X)| X € Qi). 

在 每 次 迭代 中 ， 问 题 P1(X) 的 上 界 UB 为 UB = Yo(32*)， 其 中 (2*) 表 示 当 前 已 知 的 关 
于 问题 (P1(X)) 的 最 好 可 行 点 。 

算法 的 具体 步骤 如 下 : 

步骤 1 初始 化 。 给 定 参 数 e > 0 及 向 量 参 数 o = {0,1}N， 令 k= 二 0，Qo = X? ( 即 问题 的 初 
始 向 量 空间 )， 上 界 UB = coo， 可行 点 集 F =). KE LRP (X9) 得 最 优 值 8(X?) 及 最 优 解 zo = 
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zZ(Xo)， 则 取 下 界 为 po = B(X9)。 若 z0 对 P1 是 可 行 的 ， 则 置 已 = {zx} 和 UB = Yo(z?)。 
车 UB < Bo +s， 则 算法 停止 ，z0 是 P1 的 最 优 解 ， 即 z 是 问题 LMP HRR, Vol?) W 
题 LMP 的 最 优 值 ， 否 则 令 & = 1， 转 为 执行 步骤 2; 
步骤 2 细 分 。 根 据 对 分 法 将 选 定 的 XX* 平 分 为 两 部 分 XX*t1, XFU, WX" = (XEHLXEH, 
步骤 3 可行 性 检验 。 对 每 个 子 超 矩 形 X eX, HA 


v} :二 min V; (a; X,c), 
如 果 存 在 7 e 0,1,… milf vL E UL > UB 或 中 > Inj (j 4 0). WAX 中 删除 天， 
Bi x^ =X NX} 

步骤 4 EFR URX^ 0, MERX € X" 计算 LRP (X) 88] (X) A(X) 3$ UB— 
B(X) <e WX 2X x). BW, EY 


UB := A(X), F:= FU{z(X)}, zbest .— arg min Vo(z), 


此 时 Qi := (Qi\{X*}) UX", BEEF FE B, = goin P(X); 

步骤 5 定 上 界 。 选 取 X* 的 中 点 zmid,， 车 zmid 对 P1 是 可 行 的 ， 则 已 := FU {zmid}, 4 
上 界 UB := min(Vo(z)|re F} F#0, WA 2-— arg min Volz); 

SERO KARR. WQui-QuWMX:UB-8(X) «e X EQ S Qua = $8， 则 算法 
停止 ，UB 是 问题 P1 (X0) 的 最 优 值 和 至 是 最 优 解 ， 从 而 人 是 问题 LMP (X9) 是 最 优 解 。 否 则 ， 
选取 一 个 活动 节点 XH 满足 


k+l : k+l k+1 
X mg mm B(X), X^" :-—z(X""), 
令 上 := 上 十 1 返回 步 又 2。 


4 ”算法 收敛 性 及 算 例 


如 果 算 法 不 在 有 限 步 内 结束 ， 则 很 容易 证 明 该 算法 至 少 产 生 了 一 个 无 穷 幅 套子 序 
列 {XK} (如 文献 [11j)， 即 对 所 有 的 上 有 XX*+t! c Xr*。 因 而 根据 文献 [11] 可 知 ， 存 在 x* € 
RN ext = 2*。 从 而 ， 我 们 有 如 下 重要 的 收敛 定理 。 

定理 4.1 假定 原 问题 LMP 的 全 局 解 存 在 ， 则 算法 或 者 在 有 限 步 内 求 得 LMP 的 全 局 最 优 
解 ， 或 者 算法 产生 的 无 穷 点 列 {z*} WERA z* 必 为 问题 LMP 的 全 局 最 优 解 。 

证 明 车 算法 有 限 步 终止 ， 则 结论 是 显然 的 。 因 此 假定 算法 产生 了 一 个 无 穷 迭 代 序 
列 ， 不 失 一 般 性 ， 不 妨 设 由 对 分 法 分 割 过 程 产生 一 个 无 穷 超 和 矩形 序列 X*， 对 所 有 的 上 满 
A OX*H C XK。 根据 算法 得 


< min V x* ko (XK) c XF. 
B. < min Vo(z), € erg min B(X), r'-z(X)eX 


ES (z*) 包含 于 紧 集 X09， 所 以 存在 一 个 收敛 子 序列 {23} € {rf} 设 lim z* = 人。 因而 根据 
算法 ， 存 在 一 个 下 降 超 矩形 序列 {X7} C {X} HPX EQ ER 


z EX", B,-—f(X") = WL(a"), Jim X* = (8). 


650 I E 数 学 学 R 第 27 卷 





由 定理 2.2， 有 
Jim br = lim Yor") = lim Wo(z") = Vo(2). 
因而 只 要 证 分 是 问题 LMP (X9) 的 可 行 点 。 
HFX 是 闭 的 ， 因 而 人 e X90。 下 面 利用 反 证 法 证 明 征 ;(2) < In6;, j = 1,2,… ,m。 反 设 
存在 he {1,2,… ,mj} IER U(E) > In 6h。 因 为 VL 是 连续 的 且 根 据 定理 2.2， 所 以 {VL(z"; X", 
o) 收敛 到 亚 ,( 人 2)， 必 存在 x*， 使 得 所 有 的 r+ > 地 成 立 


[Vr (z^; X",0) — V4(2)| < Vi,(£) — 1nós. 
因此 ， 对 所 有 的 > >r, AUL; X, o) > In6h， 这 上 暗示 LRP (X7) & 8] f$, 3X 5a" = 
z( X°) 矛盾。 因而 公 是 问题 LMP (X9) 的 可 行 点 。 


为 了 验证 这 个 全 局 优化 算法 ， 考 虑 下 面 两 个 LMP 问题 。 
例 4.1 图 


min (zi 十 2Z2)(Z1 — T2 + 7) 


271 + z2 X14, zy + z2 < 10, — 4zı + T2 € 0, 
st. FS = 
22; + £2 26, Tı — T2 «3, zı + 2x2 26, zı € 5. 
首先 ， 对 1; = 1, 2 求解 如 下 4 个 线性 规划 
min/max zi 
s.t. (z1, £2) € FS, 
得 到 初始 超 矩 形 Xo = {1 < zi < 5,1 < r2 < 8}， 然 后 取 o = (,0)7, e = 10-4， 用 本 文 方法 
求 得 最 优 值 为 10， 最 优 解 为 zi = 2, 23 = 8， 和 迭代 次 数 为 3 次 。 
例 4.2 
min G(z) = (z1 + 2x2 一 0.5)(2zl — z2 + 3)(z1 一 Z2 + 4.5) 
s.t. (3zl — z2 + 4)(2z1 一 1.5zz 十 3)( 一 Zi 十 2zz + 2.5) < 150, 
(一 Z1 222 十 2)(2z1 一 Z2 十 2)(1.571 4-222 一 2)( 一 Zi 十 3zs + 2.5) < 165, 
z1 2223, 
zı € [1,3], x€ [1,3]. 


Wo = (0,0)7, € = 10-4， 用 本 文 方法 求 得 最 优 值 为 166.993， 最 优 解 为 zl = 1.7473, 12 = 
1.2527， 迁 代 次 数 为 20 次， 计算 结果 说 明 该 方法 是 可 行 和 有 效 的 。 
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A Global Optimization Algorithm for Linear 
Multiplicative Programming 
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University, Changsha 410075; 2- Department of Applied Mathematics, 
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Abstract: In this paper a global optimization algorithm is proposed for the linear multiplicative pro- 
gramming (LMP). (LMP) has a number of important applications in various areas such as production- 
transportation, plant layout design, VLSI chip design. Firstly, (LMP) is transformed into an equivalent 
problem (P1). Second, utilizing the parametric linearization relaxation method, linear underestimates 
of the objective and constraint functions, linear relaxation programming (LRP) for linear multiplicative 
programming (LMP) is established over a rectangle, and the proposed deterministic global optimization 
algorithm is proposed that can converge to the globally optimal solution of (LMP). And finally the 
numerical experiment is given to illustrate the feasibility and efficiency of the proposed algorithm. 


Keywords: LMP; parametric linearization relaxation; global optimization 
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